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ABSTRACT 

We re-analyse the veloci ty-dispersion profile cr (r) at radii r > lOkpc in the Galactic stellar 
halo, recently derived bv'Battagl ia et alJ J2005h . who concluded that, for a constant velocity 
anisotropy of the tracers, these data rule out a flat circular-speed curve for the Milky Way. 
However, we demonstrate that if one makes the natural assumption that the tracer density is 
truncated at rt > 160 kpc and falls off significantly more steeply than r"^^ at r > 80kpc, then 
these data are consistent with a flat circular-speed curve and a constant velocity anisotropy 
comparable to that observed for halo stars in the Solar neighbourhood. We also consider a 
more detailed mass model with an exponential stellar disc and an extended non-isothermal 
dark-matter halo. In this two-component model, the Milky Way's virial radius and mass are 
Tvii- - 200 kpc and Mvh- ^ 1.5 x 10'^ Mq. Still assuming the tracers' velocity anisotropy to 
be constant (at y6 ^ 0.5) we again find good agreement with the observed cr(r), so long as 
the tracer density is truncated near the virial radius. These data by themselves do not allow to 
differentiate between different dark-halo or total-mass models for the Milky Way, nor between 
different velocity-anisotropy profiles for the tracers. 
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1 INTRODUCTION 

The case for dark-matter haloes around spiral galaxies is predom- 
inantly based on the flatness of the rotation curves observed for 
gas out to large radii. For large elliptical galaxies at the centres of 
clusters, similarly strong evidence comes from X-ray haloes and 
from the kinematics of stell ar tracers such as globular clusters out 
to a f ew eff"ective radii (e.g.. lC6te et alJl200lLl2003l : lRichtler et alj 
l2004h . For more isolated, 'field' ellipticals, however — or even those 
off-centre within clusters — the picture is somewhat less clear. For 
example, the projected velocity-dispersion profiles, o"(r), of plan- 
etary nebulae in fou r such galaxie s have bee n measured to sev- 
eral ef fective radii bv lMendez et alJ i200ll) and lRomanowskv et alJ 
i2003h . Although these data probe regions that are expected to be 
dynamically dominated by dark matter, the dispersion profiles are 
found to decrease outwards in a manner consistent with the as- 
sumption that light traces mass, i.e., that there are no extended 
dark haloes. However, unlike rotation curves, a-(r) profiles do not 
directly measure the total mass distribution, but are also affected 
by the density and velocity-anisotropy profiles of the tracers them- 
selves. The importance o f velocity anisotropy particularly has bee n 
emphasised by Mamon & LokasI j2005l) and iDekel et alj <2005h . 
who showed that falling dispersion profiles for elliptical galaxies 
are perfectly consistent with ma ssive dark-matter haloes, i n con- 
trast toJhe_ori^indconcliK bv lRomanowskv et alJ i2003h . 

iBattaglia et alj i2005h have recently derived a velocity- 
dispersion profile for tracers at large radii in the Galactic stellar 
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halo (globular clusters, horizontal-branch and red-giant stars, and 
dSph satellites), finding a falling (T{r) much like in the isolated el- 
lipticals. In analysing these data, they carefully consider the effects 
of velocity anisotropy and conclude that (1) a strictly flat circular- 
speed curve, corresponding to an isothermal total mass profile for 
our Galaxy, can be ruled out if the tracers' velocity anisotropy is 
spatially constant, and (2) a massive dark halo is consistent with 
the data only if its density falls off rapidly beyond r > 100 kpc 
and/or the stellar tracers are on predominantly near-circular orbits 
at these radii. 

Underlying this argument by ' Battaglia et alj is an assump- 
tion that the stellar tracers follow a power-law density profile, 
p(r) oc r^^-^ all the way out. In their Appendix B, these authors 
briefly addressed some possible consequences of relaxing this as- 
sumption, if the tracer velocity-anisotropy profile is assigned one 
specific form. However, they did not pursue any detailed mod- 
elling with the tracer p(r) allowed to deviate from a pure power 
law. The purpose of this paper is to develop such models. We find 
that, since the velocity data from IBattaglia et al] reach more than 
halfway to the virial radius of the Galaxy — probing regions where 
the stellar density is unconstrained observationally, and even ap- 
proaching what might be viewed as a natural 'edge' to the stellar 
halo — our ignorance of the detailed density distribution of the far 
stellar halo is at least as important as uncertainties in the tracers' 
velocity anisotropy when attempting to use the observed cr(r) to 
infer anything about the mass profile of the Galaxy. 

Before considering the Milky Way data specifically, let us first 
review the factors affecting the velocity dispersion of a tracer popu- 
lation. We assume dynamical equilibrium and spherical symmetry. 
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denote the tracers' density profile by p(r), and characterise the to- 
tal galactic mass distribution by its circular-speed curve V^(r) = 
GMiot{r)/r. Then the radial component iT,(r) of the tracer popula- 
tion's velocity dispersion must satisfy the Jeans equation 



2/? 



(1) 



where p(r) = 1 - crl/crj. is the usual velocity anisotropy parameter 
of the tracers (J3 = for isotropy, < /3 < 1 for radial anisotropy, 
and yS < for a tangentially biased velocity distribution). 

In the simple situation of constant anisotropy /3, and power 
laws for both the tracer density profile and the galaxy's circular- 
speed curve, p oc r"'' and Vc r", the Jeans equation requires that 
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which reduces to equation (B2) of iBattaelia et alJ i2005h in the 
special case a = (V^ = constant). Thus, both components of 
cr are proportional to Vc, and they fall with radius if and only 
if the circular speed does. More realistically however, /? and/or 
y = -dlnp/dlnr will be functions of radius; but we expect that 
at relatively large radii these will vary slowly enough for equation 
J5) to be still approximately valid and provide some heuristic in- 
sight into the general behaviour of o-(r). 

First, then, let /3(r) decrease with radius, so that the tracer 
velocity distribution becomes more tangentially biased at larger 
radii. In this case, the denominators in equation ^2) increase, so 
that (Tr(r) falls relative to V^^ir). Thus Battaglia et al. (2005), work- 
ing with strictly constant y, conclude that a rather strongly de- 
creasing anisotropy profile is nece ssary to reconcile their observed 
falling cr,.(r) profile with a standard lNavarro. Frenk&Whit^i 19961 
'NEW') model for the mass distribution of the Milky Way. 

For 7 > 2(1 + 0-) (which holds in the Galactic halo and the 
outskirts of ellipticals), the efi"ect on erg is dominated by the nu- 
merator in equation such that a decreasing fi{r) results in crg(r) 
rising relative to Vc(r). The effect on the projected line-of-sight ve- 
locity dispersion is dominated by erg (but complicated by the pro- 
jection) and results in the well-known degeneracy between the ve- 
locity anisotropy and the mass profile, which has plagued the in- 
terpretation of the velocity-dispersion profiles of ellipticals (for a 
recent exploration of this problem, see|Mamon^& Lokas 2005). 

Second, suppose instead that y(r) increases with radius (corre- 
sponding to p(r) decaying faster than a power law) while (5 remains 
constant. In this case, equation J2) tells us that both Crir) and crg{r) 
fall relative to V^ir). In fact, if — for example — the tracer density 
decays exponentially, the anisotropy is fixed at j8 = 1/2, and the 
galaxy circular speed is constant (a = 0), then we find from an 
exact solution of the Jeans equation Q that cr,. oc r"'/^: the veloc- 
ity dispersion declines like that for a power-law p(r) in a Keplerian 
potential, even though the rotation curve is flat! 

An even more drastic behaviour results from a tracer popula- 
tion which is truncated at a finite radius r,. In this case, 7 — ^ oo as 
r ^ r, and from equation J5} we expect both components of the 
velocity dispersion to vanish in this limit. Thus, the falling tr,.(r) 
observed for the stellar halo of the Milky Way may still be consis- 
tent with a roughly flat circular-speed curve, if the tracer population 
efl'ectively dies out at some large radius. This is the possibility that 
we pursue in more detail in Section|5| Note that for elliptical galax- 
ies the tracer density is known, so that this problem does not arise 
in the interpretation of their velocity dispersions. 



2 MODELLING THE MILKY WAY 
2.1 A geometrical correction factor 

The data collected by iBattaglia et alj i2005h comprise heliocen- 
tric radial velocities corrected for the Solar motion relative to the 
Galactic centre, yielding velocities, vqsr, in a Galactic standard of 
rest. However, these velocities must still be referred to the Galactic 
centre before comparing to any model. Thus, using spherical polar 
coordinates with the centre of the Milky Way at the origin, and the 
Sun on the z-axis at distance Rq = S kpc, we express the instanta- 
neous distance d from an object to the Sun as 



^r'- +Ri- 2rRo cos 



(3) 



The time derivative of this is d = vqsr = flrV, + agVg, where 

a,. = {r - Rocosd)/d and ag = (i?o sin6)/rf . (4) 

The velocity dispersion at a given Galactocentric radius r is esti- 
mated as the rms velocity averaged over all solid angles: 

c^GSR = <«?>o-? + <«e>c^e' (5) 

where (.) denotes the averaging over the sphere. Using a^. + a^ = 1 
and cTg = (1 - p)cr^, we thus obtain 

o-gsr('-) = Vl -P(r)A(r) 

with A(r) = (aj) a monotonically declining function of radius: 



A{r) = 



4r^ Sr^Ra 
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r + R, 



r-Ro 



(6) 



(7) 



The correction factor in equation j6} becomes unity for an isotropic 
velocity distribution (f} = 0), since in this case the velocity disper- 
sion is independent of direction. For non-zero yS, the correction is 
greatest for small r, where A(r) becomes largest (A = 1/2 at r = Rq 
and A = 2/3 at r = 0), whi le for r » Rq, A <K 1 and o-qsr a a,.. 

Unfortunately, Batta glia et alj|20o3) use an incorrect formula 
to convert between cr, and ctgsr [their equation (3) corresponds 
to ^6} above], which gives a correction factor larger than unity for 
j8 = 0. This significantly afl'ects the predicted ctgsr only at radii 
r < 30 kpc. However, this is where the uncertainties in the observed 
(Tgsr are smallest, resulting in substantially erroneous for any 
model fits. 



2.2 Simple models with truncated tracer populations 

Let us now employ some simple, fully analytic models of truncated 
density profiles to demonstrate that the Milky Way ctgsr data are 
consistent with a flat circular- speed curve, i.e., an isothermal total- 
mass distribution. We assume that the density of the tracer popula- 
tion, outside of some (small) core radius, is given by 



p(r) oc 



for r < r,, 
for r > r,. 



(8) 



which is plotted in Fig.Qfor y = 3.5 and n = 1, 2, or 3. At small 
radii, say r < 0.5r,, the profiles in this graph are hardly distin- 
guishable from a pure power-law p oc r"'^, which matches the den- 
sity of the Galactic stell ar halo out to r ~ 50 kpc iMorrison et alJ 
bOOdlYannv et all2OO0ft . The abruptness of the truncation towards 
larger radii depends on th e parameter n. We have used the online 
database of Ha rrisI il996h to verify that, for n > 150 kpc, a model 
with 7 = 3.5 and n = 2 provides a good description of the density 
of Galactic globular clusters at r > 3 kpc. This is thus the model we 
use for all applications to data in this paper. 
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Figure 1. Truncated density profiles of equation js) witli y = 3.5 and 
(solid), n = 2 (dotted), or n = 3 (dashed). 
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Figure 2. Radial velocity-dispersion profile for a tracer population with 
constant anisotropy f} and truncated density js) with y = 3.5 and n = 2, 
embedded in a gravitational potential with constant circular speed {a = 0). 
For a slightly falling circular-speed curve, such as predicted for dark-matter 
haloes, o",- falls off slightly faster. 



Given equation Js}, a constant anisotropy p for the tracers, 
and a power-law circular-speed curve, Vc{r) = Voirlri)" , the radial 
component of the velocity dispersion can be obtained analytically 
from the Jeans equation Q for integer n: 

^2 = K y M (_iy,-* J^_2a-2P; x] (9) 

with X = r/r, and 



1) for ^9^0. 



(10) 



r' d j - In X for ^ = 0, 

V q 

In the limit r — > r,, cr,. — > 0, as the simple considerations of Sec- 
tionQalready implied it must. (In this limit equation j9j is numer- 
ically unstable against truncation errors, but instead one may use 



-i- [ 1 + (n -H + 2(a +/?)- 211^ ) -i- + 0(6^)1 (1 1) 
n-l-ll\ 2/n + 2 J 



with e = 1 - X > 0.) 

In Fig.|2| we plot cr,.(r) for a fiat rotation curve (a = 0) and 
y = 3.5, « = 2 (appropriate to the situation in the Milky Way stel- 
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Figure 3. Comparison of truncated power-law models embedded in a flat 
circular- speed curve of Vc = 220 km s"' to the Battaglia et al. i2005) data 
for the Milky Way. Slightly better fits could be achieved with smaller trun- 
cation radius rt or a slightly falling circular-speed curve. 



lar halo), and for various /?. Evidently, the velocity dispersion de- 
clines over much the whole range of radii, even though the density 
deviates from the power-law form only for large radii. The physi- 
cal reason is that in order to achieve a truncation in density, orbits 
that reach near n must be de-populated relative to a non-truncated 
model. These orbits, however, would have contributed at smaller 
radii with rather large v,.. Note that all the models presented in 
Fig. 121 are physical in the sense that they possess a non-negative 
distribution function of the form f{E, L) = L^^^g(E). 

In Fig. |3| we compare these simple models (after applying 
the geometric correction of Section IzTl to the a-QSRir) data of 
iBattaglia et alj j2005). The tracer density again has 7 = 3.5 and 
M = 2, and we set the truncation radius to rt = 160 kpc. We as- 
sume a constant circular speed of = 220 km s"' and constant /3 
of (isotropic), 0.25, or 0.5 (radially anisotropic). All these models 
correctly reproduce the decline of cr towards larger radii. 

At small radii (r < 40 kpc), the isotropic model (solid curve) 
cannot account for the observed dispersion, but falls significantly 
below. Since at these radii 7 » 3.5 is observationally secure, we 
have from equation ^2) that cXr ~ K/VsiS = 0.53 Vc for any 
isotropic model if V^(r) is near-flat just outside the solar circle. A 
possible explanation of the observed high ctgsr for a tracer popula- 
tion with isotropic velocities requires a circular-speed curve which 
rises like V^. oc r" ' out to r ~ 30 kpc (where then = 250 km s"'). 
Alternatively and more reasonably, the data are consistent with a 
near constant Vc = 220 km s"' if the tracer velocities at r < 40 kpc 
are slightly radially biased with p = 0.25 - 0.5, comparable to 
what has been derived for halo stars in the Solar neighbourhood 
iChiba & YoshiL199&:.Gould-2004) . 



2.3 More sophisticated models 

The above analysis shows that the IBattaglia et alJ i2005l) data are 
in fact consistent with constant velocity anisotropy for the tracers 
and an isothermal total mass distribution for the Milky Way. We 
can further ask whether the same data are consistent with a more 
detailed model explicitly including both a stellar disc — which is 
well known to dominate the total mass inside the Solar circle (e.g., 
IPehnen & Binne\lll998h . and still contributes > 10% to the mass 
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_ potential: stellar disc + dark — matter halo 
tracers: r, =200 kpc. 7=3.5, n=Z, /? = 0.5 
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Figure 4. Two-component mass model for the Milky Way, and the velocity- 
dispersion profile for a tracer population in the stellar halo. The thin, solid 
(black) curve is the total circular-speed curve Vc(r) of the model (see 
eqs. 1121 and 1131 ). The bold, soHd (red) curve running through the data 
is the predicted GSR velocity-dispersion profile for a tracer population with 
a density profile truncated at = 200 kpc (see text for details) and with 
a spatially constant velocity anisotropy /? = 0.5. The thinner, dotted (red) 
curves bracketing this are the predicted itgsr for/3 = 0.5 still, but with tnm- 
cation radii rt = 170 kpc and rt = 230 kpc. The bold, dashed (blue) curve 
above these is the o"gsr that would obtain with p = 0.5 and p(r) oc r"^-^ at 
all radii (no trun cation). Data p oints are from Battaglia et aL i2005il (filled 
circles) and from lHarril i200l|) (open squares). 



at 50 kpc — and an extended dark-matter halo which is (as per cos- 
mological A'-body simulations) non-isothermal. 

We first specify an exponential stellar disc with rotation curve 



K.dC*-) = — - — = — 



where Ma 



5 .8 X 1 0'° and = 2.4 kpc according to 



iDehnen & 'Binnev^^l998^ 

Second, we take a dark-matter ha l o from the family of models 
developed bv 'Pehn en & McLa ughlin' 1*200?). Specifically, we as- 
sume that the halo follows the scaling ph/cJh r"^'''^ (consistent 
with simulations) and has an isotropic velocity distribution at its 
centre (but may be anisotropic elsewhere). Such a halo has a den- 
sity cusp pi, — > r"''' in the limit r — > 0, and ph — * r"^''' as r — > oo. 
Thus, it has a finite total mass, and its circular-speed curve is 



First, in a 'concordance' A-CDM cosmology with n„, = 0.3 and 
Qa = 0.7, the virial radius of a halo is that within which the aver- 
age density 3M(,(rvi|.)/4;rr^.^, is equal to 337pc (e.g.. Bullo ck et all 
2001). With Hq = 70km™-' Mpc"', ^ 200 kpc and Mv,,- ^ 
1.5 X 10'^ Mq for our model, both of which are reasonable for an 
L, galaxy like the Milky Way. Alternatively, haloes are commonly 
measured by the radius rim within which the average density is 
200pc. In our case, rim - 250 kpc and Miqq ^ 1.75 x 10'- Mq. For 
comparison with numbers given by Battaglia et al. ( 200^, the total 
mass inside 120 kpc is M{r < 120 kpc) = 1.05 x 10'^ Mq. 

A scale of interest in connection with numerical simulations 
of dark-matter haloes is the radius r_2 at which the local logarith- 
mic slope of the density profile is d In ph/d In r = -2. In our model, 
r_2 = (11/13)'''' ro = 27.8 kpc, an d the ratio ^200/^-2 = 9 is nicely 
consistent with the values found bv lNavarro et alJ i2004) for simu- 
lated haloes with masses in the range M200 = I - 2 x 10'^ /i^'Mq. 

To predict the kinematics of tracers in the stellar halo, we use 
the density model for the tracers that has already been employed in 
section lT^ i.e., equation (S) with 7 = 3.5 and n = 2, although we 
set the truncation radius ri = r^^ = 200 kpc, somewhat larger than 
the value used in Figure|3| We then solve the Jeans equation Q for 
an assumed spatially constant velocity anisotropy p, and apply the 
geometric correction in eq uation j6l. This yields a model crGSR(r), 
which is compared to the iBattaslia et alJ i2005l) data to compute 
X'^- The minimum^- = 9.5 (for 10 data points) is achieved with 
p - 0.5. As was mentioned above, this slight radial bias is con- 
sistent with observations of halo stars in the Solar neighbourhood 
iChiba & Yoshii 1998; Gould 2004). 

Figure B| shows our best-fit ctgsr profile against the 
iBattaglia et alj data, plus another estimate of the stellar-halo 
velocity-dispersion profile from Harris r2001). The latter is also 
constructed from velocity data for globular clusters, RR Lyrae 
stars, and dwa rf spheroidals, and so is not independent of the 
IBattaglia et alj profile; thus, we have not used it in determining 
)p- for our models. However, this alternate profile serves to con- 
firm the overall sense of the lBattaglia et alj results (and to empha- 



sise their uncertainty at the largest radii). We have also plotted in 
Figure |4] alternate models in which the tracers still have constant 
j8 = 0.5 but are truncated at r, = 170 or 230 kpc. The bold, dashed 
curve which declines only gradually towards large radius is the 
velocity-dispersion profile assuming the tracer density profile to be 
an untruncated pure power law p oc r"'^. This is the assumption 
made by Batta glia et alj i2005h in their modelling, and it clearly 
has a dramatic — even dominant — influence on the anisotropy pro- 
files they require in order to fit the observed (TQ^^ir). 



GM^ir) GMi,, 



„4/9 



4/9 



(13) 



from equation (20f) of lDehnen & McLaughUrJ ^M^. 

To fix the halo parameters, we specify the value of the total 
circular speed, = V^j + V^^^, at two Galactocentric radii: at the 
position of the Sun, r = S kpc, we require K = 220 km s"' , while at 
r = 50 kpc we adopt = 205 km s"' . This latter value corresponds 
to a total mass of Mi+Mti = 4.9x10" Mo inside r = 50 kpc, chosen 
in order to agree with the an alysis of |Kochanek ( 1996, see also 
iRohlfs & Kreitschmanij|l988l) . These two constraints then imply 
Mh,,ot = 1.10 X 10'^ Mo and ro = 40.5 kpc. 

The total mass Mh.tot here is obtained formally by integrat- 
ing the dark-matter density profile to infinity, but physically more 
meaningful is the halo mass within finite radii corresponding to 
specific over-densities relative to the critical p^ = 3Hl/SnG. 



3 DISCUSSION 

Comparing Figures|3|and|4] it is clear that in the latter we are able 
to describe the observed ctgsrC'") profile with a larger assumed trun- 
cation radius for the stellar-halo tracers than in the former, and 
also that a constant y6 = 0.5 predicts slightly lower velocity dis- 
persions at small radii in our (disc-l-halo) mass model than in the 
constant- model. These points simply reflect that the total circu- 
lar speed given by equations <12t and <13> decreases monotonically 
with radius for r > 6.6 kpc (which adds to the effect of a truncated 
tracer density in driving the decline of ( Tgsr), and in fact is l ess 
than 220 km s ~' at all radii covered by the lBattaglia et alj i2005h or 
iHarrislJioOll) data (so for fixed 7, a slightly higher y6 is required to 
give the same normalisation to the tracers' cr^, as can be seen from 
eq. Ill)- 
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Even so, all of the anisotropic models illustrated here fit the 
data equally well, and thus the tracer velocity-dispersion profile 
alone cannot be used to distinguish between different (reasonable) 
models for the total mass distribution in the Milky Way. Conversely, 
beyond concluding simply that the observed o"gsr('') profile is con- 
sistent with a tracer density profile that falls steeply at r > 80 kpc, 
we cannot use these kinematics to infer anything more detailed 
about the behaviour of p{r) for the far stellar halo. 

Nor can strong constraints be placed on the velocity 
anisotropy of the tracer population. In the context of a spatially 
constant j8, as we have assumed, all that can be said with any con- 
fidence is that a slight radial bias (Ji ~ 0.2 - 0.6, depending on 
the Galactic mass model) is required to explain the amplitude of 
(Tgsr at r < 40 kpc. If the assumed form of the total mass dis- 
tribution deviates significantly from the models with constant or 
slowly- varying V^ir) that we have explored, then a variety of p{r) 
behaviour is likely allowed by the data. 

Clearly, untangling the degeneracies between the total M(r) 
and tracer p(r) and j8(r), to put quantitative limits on any one pro- 
file from knowledge of the others, is a daunting task that will re- 
quire much closer attention to a wider variety of data beyond just 
o"gsr('')- For example, the models for V^ir) that we have worked 
with here are extremely simple and would surely require some 
modification in detail if we attempted to take accurately into ac- 
count the many other constraints on the Galactic mass distribution. 
But more detailed modelling will also have to allow for the fact 
that th e tracers contributing to the Battaglia et al. 1 2005) and Harris! 
i200lh velocity-dispersion profiles are very much a 'mixed bag': 
the halo globular clusters, for instance, may not have the same 
anisotropy profile as the field RR Lyrae and red giants; and the 
satellite dwarf spheroidals likely do not follow the same density 
profile as the stars and globulars. Ideally, we need information on 
separate crir), p(r), and fi(r) profiles for each of the different tracer 
populations. Current data simply do not provide this. 



4 SUMMARY 

We have demonstrated in this study that the falling velocity dis- 
persion found by Battaglia et al. (2005) in the outer stellar halo of 
the Milky Way is consistent with a constant (and reasonable) ve- 
locity anisotropy for the stellar tracers, and either a perfectly flat 
circular-speed curve (i.e. an isothermal total mass distribution) or 
a standard CDM dark-matter halo combined with an exponential 
stellar disc. By contrast, Battaglia et al. argued that their data were 
inconsistent with and /? both being strictly constant, and that they 
required a strongly varying (and r ather unusual) /?(?•) pr ofile to be 
made compatible with the common iNavarro et alljl99d) model for 
a dark-matter halo. ' 

As we have emphasised, the primary reason for the difference 
between our conclusions and theirs is that we allow for the tracer 
density to die out at r > 160kp c, close to the Milky Way's virial 
radius, whereas Bat taglia et al J assume that p oc r^^-^, which is valid 
for r < 50 kpc, continues to hold for all larger radii. An incorrect 
tr ansformation fro m modelled (Tr{r) to observed crGSR('') on the part 
of lBattaglia et alj (see Section im has certainly also contributed to 

' When lBattaglia et alj|20o3) discuss 'NFW' or 'truncated flat' or 'isother- 
mal' haloes, they really model the total Galactic mass distribution with 
those functions. But as we discussed in connection with Figure|3 the stellar 
disc contributes significantly to the total circular speed out to tens of kpc, 
and the total VcC") differs from that of the dark halo until ver)' large radii. 



the discrepancy, insofar as it affected the;^'^ values of their detailed 
model fits. 

As we discussed, these data by themselves are insufficient for 
differention between different total-mass models for the Galaxy, 
nor between different velocity-anisotropy profiles for the tracers. 
This is essentially because the measured cr(r) depends on both of 
the above (the well-known degeneracy between velocity anisotropy 
and mass profile) and on the tracer's density profile, which is little 
constrained observationally at r > 50 kpc. 
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